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ABSTRACT 

We discuss the results of a multiwavelength differential imaging lab exper- 
iment with the High Contrast Imaging Testbed (HCIT) at the Jet Propulsion 
Laboratory. The HCIT combines a Lyot coronagraph with a Xinetics deformable 
mirror in a vacuum environment to simulate a space telescope in order to test 
technologies and algorithms for a future exoplanet coronagraph mission. At 
present, ground based telescopes have achieved significant attenuation of speckle 
noise using the technique of spectral differential imaging (SDI). We test whether 
ground-based SDI can be generalized to a non-simultaneous spectral differential 
imaging technique (NSDI) for a space mission. In our lab experiment, a series of 
5 filter images centered around the 02(A) absorption feature at 0.762 /im were 
acquired at nominal contrast values of 10~^, 10~^, 10~^, and 10"^. Outside the 
dark hole, single differences of images improve contrast by a factor of ~6. Inside 
the dark hole, we found significant speckle chromatism as a function of wave- 
length offset from the nulling wavelength, leading to a contrast degradation by a 
factor of 7.2 across the entire ~80 nm bandwidth. This effect likely stems from 
the chromatic behavior of the current occulter. New, less chromatic occulters 
are currently in development; we expect that these new occulters will resolve the 
speckle chromatism issue. 



1. Introduction 



While over 300 exoplanets have now been discovered indirectly via radial velocity, tran- 
sits, or gravitational lensing, to date only a handful of planets have been directly imaged 
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around a main sequence or pre- main sequence star (IKalas et al.ll2008l : iMarois et al.ll2008l ). 
Most other detection methods (e.g. radial velocity and transit) are more sensitive to close-in 
planets than those at greater distance from their stars. Direct imaging is the only method 
currently capable of detecting planets >10 AU from their parent stars. Numerous adaptive 
optics and space based surveys to directly image young Jup iter mass planets in the near-IR 
( iBiller et al.l 120071 : iLafreniere et al.l 120071 : iFarihi et al.l 120071 ) have recently been completed; 
these surveys have placed strong statistical limits on extrasolar planet distributions but have 
produced no planet detections. Clearly, directly imaging planets is a nontrivial undertaking 
- and one that becomes even more difficult for terrestrial planets. 

In order to successfully image a planet one must first overcome: 1) the huge contrast 
difference between star and planet (the very youngest Jovian planets will be 10~^ times 
fainter than their parent star, while a terrestrial planet will be ~10^^'^ times fainter). Mod- 
ern approa ches to planet-finding of ten employ a Lyot coronagraph to achieve these con- 
trasts (e.g. iKuchner fc Traubll2002l ). 2) Even once the initial star- planet contrast has been 
achieved, residual superspeckles remain (from errors in the wavefront due to imperfect op- 
tics) and produce a limiting "speckle" noise fioor which prevents planet detection. This is 
seen both from the ground and space (e.g. the "breathing" effect seen by HST). 

A number of instrumental speckle suppression techniques have been suggested and im- 
plemented on ground base telescopes (and in testbed scenarios for space base d telescopes), 



including azimuthal differe ntial imaging (aka ADI or roll subtraction, see e.g. iMarois et al. 
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Lenzen et al.ll2004 : Chun et al.ll2008l) . imaging polarimetry ( Perrin et al 
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Qppenheimer et al.ll200£ ), var i ous methods of focal pl ane wave 



20081 : iKenworthv et al.ll2006f ). and the 



use of Integrated Field Spectrographs (IFS) t o separate speckles (which modulate with wave- 



length) from real c ompanions (which do not) (ILarkin et al.ll2006l : lMcElwain et al.ll2007l . 12008 



Claudi et al.ll2008l ). In this paper we focus on the possible application of the SDI technique 
to a space-based telescope platform. 

From the ground, the SDI technique has been used after adaptive optics (AO) correction 
to remove residual speckle noise witho ut removing planet light. SDI has been implemented 
on the TRIDENT device at the CFHT (Marois et al. 12005), on the Simultaneous Differential 



Imager at the VLT and MMT (IBiller et al.ll2007l : iLenzen et al.ll2004j ) and, more re cently, with 



a cor onagraph, on the Near-Infrared Coronagraphic Imager at Gemini South (jChun et al. 



20081 ). SDI attenuates speckle noise by exploiting a spectral feature in the desired target 
(previous implementations utilize the 1.6 iJ,m methane absorption feature robustly observed 
in substellar objects cooler than 1400 K). Images are taken simultaneously both within and 
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outside the chosen absorption feature. Due to the simultaneity of the observations, the star 
and the coherent speckle pattern are largely identical in both filters, while a faint companion 
with that absorption feature is bright in the continuum filter and faint in the absorbed filter. 
Subtracting the absorbed filter image from the continuum filter image thus removes the 
starlight and speckle patterns while a real companion with the chosen absorption feature 
remains in the subtraction image. In other words, the off-absorption band image acts as 
an ideal reference point spread function (henceforth PSF) for the absorption band image. 
Utilizing a signature spectral feature of substellar objects greatly reduces the number of 
false positives detected (e.g. a faint background object, while real, will drop out of the 
SDI subtraction since it will not have methane absorption). From the ground, significant 
attenuation of speckle noise has been achieved with SDI. For example, iBiller et al.l (120071 ) 
achieved H band contrasts >25000 (5 a AF1(1.575 yum) > 10 mag, AH > 10.8 mag for a T6 
spectral type) at a separation of 0.5" from the primary star. With this degree of attenuation, 
it is possible t o image (5 a detection) a 2 Jupiter m ass planet at 5 AU around a 10 Myr MO 
star at 10 pc jBiller et allboO?! : iFortnev et aPboOsh . 



SDI on ground based telescopes provides significant speckle attenuations down to star- 
planet contrasts of ~0.3-lxl0~^. To test the classical SDI technique at contrasts of 10~^ 
to 10~^ and to determine whether such a technique would be applicable for the detection 
of terrestrial planets for a possible exoplanet coronagraph, e.g. a Terrestrial Planet Finder 
mission, we implemented a similar multiwavelength differential imaging scheme for the High 
Contr ast Imaging Testbed (h enceforth HCIT, described in detail in the supplemental mate- 
rial to Trauger fc Traubl[2007h at the Jet Propulsion Laboratory. The HCIT combines a band- 
limited Lyot coronagraph with a 32x32 actuator Xinetics deformable mirror in a vacuum 
environment (ITrauger et al.ll2004j ) to simulate a space telescope in order to test technologies 
and algorithms for future coronagraph missions. The HCIT has demonstrated extremely high 
contrasts and very fine wavefront control ; contrasts of 10~^° hav e been achieved, sufficient 
to image an Earth-like extrasolar planet (ITrauger fc Traubl 120071 ). 



2. Experimental Design 



Five 2% bandwidth filters were selected near the prominent 02(A) absorption feature at 
0.762 /zm (seen in Eart h's atmosphe r e and expected for any terrestrial extrasolar planet with 
an oxygen atmosphere. IWoolf et al.l (120021 )) with approximate central wavelengths of Fl(768 
nm), F2(784 nm), F3(800 nm), F4(816 nm), and F5(832 nm). Exact filter wavelengths and 
bandwidths are presented in Table [1] - however, approximate wavelengths are used to refer 
to the filters throughout the text. Two sets of images were taken in each filter - a long 
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set, with 90 s exposure time per filter, and a short set, with 5 s exposure time per filter. 
The filter set was spaced across a considerable bandwidth in wavelength in order to measure 
speckle chromatism as a function of wavelength. The light source used for this experiment 
is a passively Q-switched Nd:YAG 1060 nm laser. This laser delivers subnanosecond pulses 
into a single-mode supercontinuum photonic crystal fiber, producing light with a relatively 
smooth spectrum from 600 nm to 1600 nm and prov iding sufficient bandwidth and stability 
to conduct this experiment (ITrauger fc Traubl 120071 ). 



For ground based observing with adaptive optics (AO), quasi-static superspeckles due to 
residual instrumental errors remain even after AO correction ( Marois et al.lEoOS : Masciadri et al 



20051 ). These superspeckles are stable on timescales o f minutes to hours, but still vary quickly 



enough over a typical observation (1-2 hours, see e.g. iHinkley et al.ll2007l : iMarois et al. 



2006f ) 



that a reference PSF must be built on similar timescales in order to overcome the stochastic 
speckle noise floor. (To build this PSF, ADI uses rotation on the sky to decorrelate speck- 
les from real object, SDI uses wavelength diversity, and imaging polarimetry uses images 
taken in different polarization states.) Thus, on the ground, simultaneous imaging in sev- 
eral bandwidths is desirable to successfully implement spectral differential imaging (SDI). 
For space-based observing, however, speckles are stable on m uch longer timescales of hours 



to days, making simultaneity of imaging less necessary (e.g. iSparks &: FordI |2002| ) . While 
proposed instrument designs (such as those for the Terrestrial Planet Finder and the Eclipse 
mission) utilize simultaneous imaging via either an integral field unit or dichroic beamslitter, 
no multispectral backend imager was available at the HCIT. 



Trauger fc Traubl (120071 ) demonstrated extreme speckle stability at the HCIT during 
a "movie experiment" where 480 individual snapshot images were recorded over a 5 hour 
period. They found that the change in speckle contrast in HCIT during this period (where 
the apparatus was allowed to drift freely) was ~0.1xlO^^° per 5 hours. This demonstrated 
speckle stability allowed the current experiment to be conducted non-simultaneously with 
valid results. Additionally, the use of the same imaging channel avoids non-common path 
wavefront errors between images and thus allows for a cleaner determination of inherent 
speckle chromaticity. The present multi-wavelength differential imaging experiment mea- 
sures speckle evolution as a function of wavelength and contrast level. We test whether the 
ground-based simultaneous spectral differential imaging technique can be generalized to a 
non- simultaneous spectral differential imaging technique for a space mission. 

By using a 5 filter set, we can attempt to correct for speckle chromatis m using both 



single differences of images and the double difference technique pioneered by iMarois et al. 



( 120001 ). Since the radial position of the speckle pattern in each filter image is proportional to 



, the platescale of each image must be scaled so that the speckles in each filter fall at the 
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same radii despite chromatic differences. After this scahng, a number of single differences as 
well as a double difference of images in the Fl through F5 filters are calculated: 

/31 = (F3 - Fl) (1) 



/32 = (F3 - F2) 



(2) 



J34 = (F3 - F4) 



(3) 



/35 = (F3 - F5) 



(4) 



/4332 = (F4 - F3) - (F3 - F2) 



(5) 



'5331 



(F5 - F3) - (F3 - Fl) 



(6) 



In this document, we discuss preliminary results using the classical SDI data reduction 
method (i.e. no Fresnel propagation is considered.) 



3. Data Acquisition and Reduction 

Data were acquired in January 2007. A series of 5 filter images were acquired at each 
nominal contrast value (10~^, 10~^, 10~^, and 10~^). High contrasts were achieved via 
speckle nulling in the F3 800 nm filter. Speckle nulling iteratively removes speckles locally 
by correcting the wavefront with the deformable mirror. Speckles are "dialed o ut" one by one 



by programming in the appropriate antispeckle with the deformable mirror (IMalbet et al. 



I995I : iTraueer et al.ll2002l : iBorde fc Traubl I2OO6I : iTraueer et al.ll2006f ). The HCIT operates 
at such high contrasts that speckles due to both phase errors and amplitude errors in the 
wavefront must be corrected. (From the ground, at contrasts of ~10~^ - 10"^, p hase er- 
rors dominate so cornplete ly that amplitude errors can effectively be ignored, see e.g. iGuyon 
2OO5I : ISoummer et al.l 120071 ). However, only phase speckles can be dialed into the deformable 
mirro r. This creates a symmetry issue when removing amplitude speckles ( IBorde &: Traub 
20061 ). In the image, both phase and amplitude speckles will appear as a pair of symmetric 
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point sources around the central PSF. However, when inspected in image plane electric field 
instead of image amplitude, phase and amplitude speckles appear rather different. Phase 
speckles are antisymmetric, while amplitude speckles are symmetric. Adding the opposite 
phase speckle on the DM will completely cancel out a phase speckle, but will only cancel out 
one side of an amplitu de speckle. The other side of the amplitude speckle will be intensified 



( jBorde fc Traubll2006l ). In our optical system, phase speckles d ominate to contrasts of 10 



up to these contrasts, a symmetric dark hole can be generated (ITrauger et al.ll2004j ). At con- 
trasts better than 10"'', speckles are produced by a mix of phase and amplitude aberrations 
in the wavefront and the dark hole becomes asymmetric. 

In this manner, the scattered light producing the speckles can be removed locally, how- 
ever, the speckle solution is inherently asymmetric; high contrasts are achieved only in the 
right-side dark hole region. (Fig. [1] shows a schematic of the darkest portion of the dark hole 
region on the right and a comparison region on the left.) This is quite clear in our images 
(Fig. [2] to Fig. [5]) - at a contrast of 10~^, a slight dark hole is seen on both sides of the 
image but at contrasts of 10"'' or better when amplitude errors must be taken into account, 
the dark hole appears only on the right side of the image. The size of this dark region is 
determined by the control radius of the deformable mirror, which is in turn determined by 
the total number of DM actuators and actuator spacing. 

Once each contrast level was achieved via speckle nulling, two images (with exposure 
times of 90 s and 5 s) were acquired at each of the five filter wavelengths. Per-pixel S/N 
ratios (for pixels representative of the average contrast value in the image dark hole) are 
presented in Tab. [2l The gain of the CCD was 2.5 electrons/DN, with a read noise of 6 
electrons/pix rms. The total noise per pixel is a combination of shot noise and read noise: 



N = ^jNl,, + Nl^, (7) 

We focus on the 90 s dataset, since its signal to noise ratio is considerably higher (by a 
factor of >4) than that of the 5 s dataset. Unfortunately, significant saturation within the 
dark hole region (~5% of pixels) occurred in the 10~^ contrast images. We display these 
images since their qualitative features are still of interest but exclude them from further 
analysis. Some saturation is seen in the 10"'^ to 10~^ contrast images, but only far outside 
the dark hole and comparison regions. In fact, for these images, only ghost point sources 
considerably outside the dark hole were saturated. (These ghost point sources are produced 
by aliasing outside the control radius of the DM and are 31 A/D from the occulted star 
- in other words, far from any analysis region). Coronagraph images were bias and dark 
corrected, but no flat flelding was p erformed, as the CCD already exhibits a pixel-to-pixel 



DQE uniformity of better than 1% (ITrauger &: Traubl 120071 ). Since the images in the nulled 
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wavelength are of the highest quahty and are the most rehable, the F3 800 nm filter image 
was used as the "master" image (for adjustment of platescale and ahgnment purposes) in 
this analysis. 

Fluxes in each wavelength image were converted to the equivalent contrast via the 
following steps: 1) a "stellar point spread function" (henceforth PSF) was derived by by 
offsetting the coronagraph mask by 144 fim to the first maximum in transmission (in other 
words, moving the star off the mask) and taking a 30 /xs snapshot image, 2) the "stellar 
peak" was estimated from the snapshot image and scaled to the appropriate peak signal 
for a 90 s data image, 3) each 90 s image was divided by the estimated peak signal to 
convert to contrast units, and 4) each peak signal corrected image was then divided by the 
analytic attenuation profile of the coronagraph mask (known to closely match the measured 
attenuation profile) to corr ect for attenuatio n . For more details on this process, see the 



supplemental material from lTrauger fc Traubl (120071 ). 



Since the radial position of the speckle pattern in each filter image is proportional to 
A, the platescale of each image was scaled to the 800 nm image so that the speckles in each 
filter fall at the same radii despite chromatic differences. The exact filter central wavelengths 
presented in Table 1 were used for this scaling. To estimate the resampling noise introduced 
by wavelength scaling, we rescaled 20 randomly chosen 40x60 pixel subimages from within 
each of the dark holes of the 10~^, 10~^, and 10~^ images (4 wavelengths x 3 contrast levels, 
so 240 subimages total were considered in this analysis). We then calculated the standard 
deviation in each subimage (using robust_sigma in IDL) before and after scaling (aorig and 
scaled respectively). Resampling noise was estimated to be: 



■^resamp StddcV {jJ Qj-ig (^scaled) (8) 

We found values for the resampling noise of 3.5-6.6x10^^^, which, even at our best 
contrasts of 10~^, are negligible. 

Finally, each wavelength image was aligned to the 800 nm image using a custom shif t- 



and-subtract alignment algorithm (alignments are to 0.1 pixel precision. iBiller et al.l (120071 )). 
The resulting aligned images were subtracted from the master 800 nm image. Galleries of 
each aligned image and the resulting subtraction are presented for each nominal contrast level 
in Fig. [2] to Fig. [51 These images are shown with a logarithmic stretch from contrast levels of 
to 10~^. Speckle rms measured in several regions outside the dark hole were reduced by a 
factor of 5-50 after subtraction (for details on speckle rms/contrast calculations, see Section 
4 below). Two double differences were also calculated and are shown for 10~^ contrast in 
the dark hole (same logarithmic stretch, also by nominal contrast level) in Fig. [61 



- 8 - 



4. Analysis 



To quantify the level of speckle attenuation available through the differential imaging 
technique both inside and outside the dark hole, we calculated the contrast levels before 
and after subtraction in two 20x60 pixel regions of the chip for all our reduced images. 
Contrast (or more simply, residual spe ckle rms) is e stima ted as the standard deviation over 
the region, an approach adopted from iBiller et al.l (120071 ) . One major caveat is attached to 
this approach - it assumes Gaussian statistics, whereas speck l es have been shown to follow a 
Rician probability density function (lAime fc Soummeii l2004l : iFitzgerald et al.l 120071 ) . Using 
Gaussian statistics w ill seriously overestimate the confidence level of a point source detection 
( iMarois et al.l l2008bl ) . However, this app r oach has been used r o bustly as a diagnos tic of 
contrast performance (ITrauger et al.ll2004j : iLafreniere et al.ll2007l : iHinkley et al.ll2007l ). For 
this testbed experiment, where we are simply concerned with speckle residuals and not 
detections of actual point sources, this method is appropriate for our use. 

Regions used for this analysis are shown in Fig. [H The dark hole region lies at 4-9 A/D 
from the occulted stellar image while the left-side "comparison region" lies 19-24 A/D to 
the left of the occulted stellar image (A/D=4.67 pixels at 785 nm for the HCIT, Trauger & 
Traub 2007 supplemental material). One major caveat exists regarding our choice of left- 
side "comparison region" - we chose a comparison region outside the control radius of the 
deformable mirror, since the speckle nulling algorithm adds speckle noise to the left side of 
the image when removing an amplitude speckle from the right side of the image. Before 
nulling we note that a radial trend in speckle intensity exists - there are more and brighter 
speckles closer to the center of each image. Thus, the speckle noise within the left-side region 
is somewhat lower than the pre-nuU speckle noise in the dark-hole region. For this reason, 
it is not instructive to compare speckle noise between the right and left side regions, but 
instead to compare the degree of attenuation provided by the differential imaging technique 
in both of these regions. Contrast information is presented in Table 2 and is plotted for the 
dark hole region as a function of AA from the nulling wavelength (800 nm) in Fig. [71 Outside 
the dark hole, the single differences improve contrast by a factor of ~6, meaning that the 
NSDI method will be highly applicable in moderate (10^^ - 10^^) contrast systems without 
a DM to iteratively remove speckles (for instance, this technique might be very appropriate 
to utilize with a coronagraphic Jovian planet imager without a DM). Inside the dark hole, 
contrast appears to depend strongly on AA from the nulling wavelength, especially at higher 
contrasts. For lower contrast levels, a similar speckle suppression with single subtraction 
is also observed within the dark hole. At high contrasts (10~^,10~^), however, apparent 
chromatic variation becomes very important and the speckle pattern appears to decorrelate. 
For instance, significantly more speckle noise appears in the 832 nm vs. the 800 nm image, 
meaning that a subtraction of these two images degrades the achieved contrast relative to 
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the 800 nm image alone. For the same reason, the double differences do not provide an 
advantage over the single differences - contrast inside and outside of the dark hole is not 
decreased compared to the single differences. 

Increased speckle noise in the dark hole as image wavelength differs from the nulling 
wavelength translates to a lower contrast in that image relative to the image at the nulled 
wavelength. Our multifilter experiment lets us simulate the variation in achieved contrast as 
a function of wavelength that we would expect within a wideband image. Our filters cover a 
wavelength range of ~80 nm - equivalent to a 10% bandwidth filter. At an 800 nm contrast 
of 10~^, contrast is degraded by a factor of 1.7 at the red end (830.9 nm) compared to the 
nulled wavelength (798.9 nm) and contrast is degraded by a factor of 2 at the blue end (768.1 
nm) compared to the nulled wavelength. Over the entire ~80 nm bandwidth, contrast is 
on average degraded by a factor of 1.5 compared to contrast at the nulled wavelength (and, 
hence, the contrast would be degraded by a factor of ~1.5x over a 10% bandwidth filter than 
over a 2% bandwidth filter.) At the highest contrasts (10~^) in the 800 nm image, contrast 
over the entire bandpass is considerably diminished - at the red end (830.9 nm), contrast is 
degraded by a factor of 11 compared to the nulled wavelength (798.9 nm) and at the blue end 
(768.1 nm), contrast is degraded by a factor of 15 compared to the nulled wavelength. Over 
the entire ~80 nm bandwidth, contrast is degraded on average by a factor of 7.2 compared 
to the nulled wavelength. Thus, apparent speckle chromatism in these images is sufficient 
to predict considerable contrast degradation in a wide band (~80 nm) filter. 

To determine the source of speckle chromatism in these images, we must take into 
account various sources of wavefront error in the HCIT. There are three main categories of 
wavefront errors we must consider: 

1) Wavelength scaling phase errors in the pupil plane. These speckles scale 
according to A and thus can be aligned across different wavelengths. Additionally, they are 
antisymmetric and can be completely removed by the deformable mirror. 

2) Wavelength scaling amplitude errors in the pupil plane. These speckles scale 
according to They can be removed on one side of the image using the deformable mirror, 
but will produce an increase in the corresponding speckle on the other side of the image. 

3) Wavelength dependent phase errors in the focal plane due to the occulter. 

In other words, each wavelength sees different occulter properties. The current occulter is 
fabricated from High Energy Beam Sensitive (HEBS) glass, which varies in transmission 
according to wavelength and thus will introduce slightly different phase errors at different 
wavelengths. (The relationship between wav elength and phase shift for HEBS glass is plot- 



ted in Figure 2 of iMoody &: Traugerl 120071 ) We hypothesize that the HEBS glass in the 
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occulter is responsible for the speckle chromatism observed in this experiment. At lower 
contrasts the phase errors due to wavelength dependent transmission in the HEBS glass are 
negligible, but at 10~^ contrasts, they become quite important. These speckles are thus 
chromatic - they do not scale by A and are likely the culprit in causing the characteristic 
"parabola" pattern in contrast seen in Fig. [71 

In Fig. [HI we compare our results to in-house theoretical predictions for this experiment 
(at slightly better contrasts of 10~^'^ and with somewhat different comparison regions) and 
find the same qualitative result. The testbed optics are modeled as a series of optical surfaces . 



each of which has an interferometrically measured surface map (see, e.g. [Fang et al.ll2005l ) 
The occulting mask is ass umed to have a profile of intensity transmission that follows a 
fourth-order sinc^ function flKuchner &: Traubll2002l ). an d has a transmitted pha se that varies 
as a function of intensity transmission, as measured by iHalverson et al.l (l2005l ). Wavefront 
correction is perform ed using a model of the speckle nulling algorithm, as described in 
Trauger et al.l (120041 ). To simulate the image recorded on the camera through an individual 
bandpass filter, a number of monochromatic wavelengths (typically at least 3 wavelengths per 
filter) are propagated through the model and summed by intensity at the camera. The model 
iterates the speckle nulling algorithm until the contrast improvements in s uccessive iterations 
become ne gligible. More details rega r ding modeling can be found in [Moody fc Trauger 



( |2007[ ) and [Moody. Gordon. &: Traugerl (120081 ). Models using the current generation of HEBS 
glass occulters can qualitatively reproduce the speckle chromaticity found in the current 
experiment. 



Moody fc Traugerl (|2007[ ) used this same modeling approach to simulate the contrast 



performance of next generation hybrid coronagraph mask designs (metallic and dielectric 
thin films deposited on a glass substrate) using a similar band-limited coronagraph setup 
and predict considerably less chromatic coronagraph performance from these masks. Thus 
we suggest that the speckle chromaticity seen in this experiment is due to the A-dependent 
behavior of the HEBS glass which ma kes up the occulter and not from th e intrinsic design of 
the band-limited coronagraphic mask. [Moody. Gordon, fc Trauger[ (120081 ) have implemented 
a metallic thin film mask at the HCIT (and achieved the best published contrasts to date at 
the HCIT with it). Hybrid masks are just becoming available at the HCIT and have already 
yielded considerable improvements in contrast over a 20% bandwidth filter (Moody, private 
communication). When these hybrid masks have been fully implemented at the HCIT, this 
experiment will be repeated to test this hypothesis. 
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5. Discussion, Conclusions, and Future Work 

Using single wavelength speckle nulling with a band-limited Lyot coronagraph and op- 
erating at nominal contrasts of 10~^ to 10~^, single differences of filter images can reduce 
speckle noise outside of the dark hole by factors of 5 - 50. For contrasts of 10~^ to 10~®, a 
similar result is also found within the dark hole, with speckle attenuation achieved through 
a single difference decreasing as a function of increasing contrast. However, at high contrasts 
(10^^), considerable increase in speckle RMS between filters is observed to "pollute" the 
dark hole in all of our single differences (where, in each difference, only one of the two wave- 
lengths had undergone optical speckle nulling). At all contrast levels, a double difference of 
images does not seem to decrease speckle noise relative to the single differences. Significant 
differences in contrast (~2-10xRMS) are found between filters separated in wavelength by 
AA >20 nm. Contrast degradation between filters increases strongly as a function of in- 
creasing contrast at the nulling wavelength. This chromatic speckle noise is likely due to the 
chromaticity of the materials of the occulter itself. 

However, this experiment is only one possible implementation using specific choices of 
coronagraphic mask/occulter and nulling algorithms available in January 2007. It is clear 
that the occulter used in January 2007 possessed significant wavelength dependent phase 
errors, causing speckle decorrelation between filters. This decorrelation was exacerbated by 
using a speckle nulling algorithm which calculated a speckle solution at only one wavelength. 
Since 2007, significant effort has gone into developing non-chromatic coronagraphic focal 
plane occulters and nulling algorithms for the HCIT. Thus, another implementation of this 
experiment using these improvements would likely be able to extend NSDI performance to 
even better contrasts. Specific improvements possible include: 

1) The major limiting factor for this experiment is likely chromatism due to the materials 
of the occulter itself. However , since January 2007, a variety of less chromatic occulter 



options are under development. iMoody fc Traugen (120071 ) explore hybrid coronagraph mask 



designs (metallic and dielectric thin films deposited on a glass substrate) which promise 
considerably less chromatic coronagraph performance. These masks are still in development 
but will in principle allow the extension of differential imaging techniques to 10^^ - 10^^° 
contrasts. 

2) Algorithms which produce a multiwa yelength speckle so lution are now available 



(e.g. multiwavelength electric field conjugation. ICive'on et al.ll2007l ) instead of speckle nulling 



at a specific wavelength. Utilizing one of these newer algorithms will produce a broadband 
solution without degradation at either end of the bandpass. 

A future iteration of this experiment will be conducted once these improvements are 
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implemented. 

Part of the research described in this paper was carried out at the Jet Propulsion Lab- 
oratory, California Institute of Technology, under a contract with the National Aeronautics 
and Space Administration. Support for B.B. was provided by NASA through the NASA 
GSRP grant NNG04GN95H and the Hubble Fellowship grant HST-HF-01204.01-A awarded 
by the Space Telescope Science Institute, which is operated by the Association of Universities 
for Research in Astronomy, Inc., for NASA, under contract NAS 5-26555. 
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Table 1. Filter Wavelengths and Bandwidths 



Filter Center Wavelength FWHM 



Fl 768.1 nm 15.1 nm 

F2 782.5 nm 15.6 nm 

F3 798.9 nm 15.9 nm 

F4 814.8 nm 16.3 nm 

F5 830.9 mil 15.5 nm 



Table 2. Per- Pixel S/N Ratio in the Dark Hole 





Contrast 




Exp. Time 


Fl 


F2 


F3 


F4 


F5 


10- 


pixel in 10"^ 


image 


90 s 


179.9 


182.4 


179.9 


187.9 


142.2 


10- 


pixel in 10~^ 


image 


90 s 


45.6 


47.1 


45.8 


68.3 


68.8 


10- 


pixel in 10~^ 


image 


90 s 


21.2 


21.5 


19.2 


19.2 


18.1 


10- 


pixel in 10~^ 


image 


90 s 


4.4 


2.5 


2.7 


2.5 


2.7 


10- 


-6 pixel in 10"^ 


image 


5 s 


42.0 


42.6 


42.0 


43.9 


33.0 


10- 


pixel in 10"'' 


image 


5 s 


9.5 


9.9 


9.5 


15.2 


15.3 


10- 


-® pixel in 10"^ 


image 


5 s 


3.4 


3.5 


2.9 


2.9 


2.7 


10- 


pixel in 10~^ 


image 


5 s 


0.3 


0.2 


0.2 


0.2 


0.2 
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Table 3. Speckle RMS in Right-Side Dark Hole (dh) and Left-Side Comparison Region 

(oh) 



Image 


Region 


10-^ 




io-« 




10-9 




Fl 768.1 nm 


dh 


7.1x10' 


Q 

-o 


1.4x10" 


o 

-o 


1.2x10" 


Q 




oh 


4.0x10" 


-8 


4.7x10" 


-8 


5.4x10" 


-8 


F2 782.5 nm 


dh 


7.0x10" 


Q 


1.1x10" 


o 

-o 


3.7x10" 


Q 




oh 


3.7x10 


-8 


4.3x10 


-8 


5.0x10 


-8 


F3 798.9 nm (nulled A) 


dh 


6.9x10" 


Q 

-o 


7.3x10" 


Q 


7.9x10- 


1 n 
■iU 




oh 


3.7x10 


-8 


4.4x10' 


-8 


4.9x10" 


-8 


F4 814.8 nm 


dh 


7.7x10" 


Q 
O 


1.0x10" 


Q 
C5 


3.0x10" 






oh 


3.6x10" 


-8 


4.2x10" 


-8 


4.7x10" 


-8 


F5 830.9 nm 


dh 


7.8x10" 


-8 


1.3x10" 


-8 


8.9x10" 


-9 




oh 


3.5x10" 


-8 


4.1x10" 


-8 


4.7x10" 


-8 


798.9 nm - 768.1 nm 


dh 


1.7x10" 


-8 


1.0x10" 


-8 


1.2x10" 


-8 




oh 


6.4x10" 


-9 


6.8x10" 


'9 


9.3x10" 


-9 


798.9 nm - 782.5 nm 


dh 


7.5x10" 


-9 


3.9x10" 


-9 


3.0x10" 


-9 




oh 


2.6x10" 


-9 


2.9x10" 


-9 


3.5x10" 


-9 


798.9 nm - 814.8 nm 


dh 


7.6x10" 


-9 


3.8x10" 


-9 


2.5x10" 


-9 




oh 


2.7x10" 


-9 


3.1x10" 


-0 


1.0x10" 


-9 


798.9 nm - 830.9 nm 


dh 


1.3x10" 


-8 


8.6x10" 


-9 


8.5x10" 


-9 




oh 


4.7x10" 


-9 


5.3x10" 


-9 


5.8x10" 


-9 


(814.8 nm - 798.9 nm) 


dh 


8.0x10" 


-9 


3.7x10" 


-9 


4.8x10" 


-9 


- (798.9 nm - 782.5 nm) 


oh 


5.2x10" 


-9 


5.2x10" 


-9 


5.4x10" 


-9 


(830.9 nm - 798.9 nm) 


dh 


1.8x10" 


-8 


1.3x10" 


-8 


1.9x10" 


-8 


- (798.9 nm - 768.1 nm) 


oh 


9.0x10" 


-9 


1.1x10" 


-8 


1.4x10" 


-8 



-17- 




Fig. 1. — Boxes used for speckle RMS calculation (Table 2) 
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Fig. 2. — Left:RAW, Right:NSDI. Gallery of single wavelength images and single differences 
with a nominal contrast level of 10~®. All images are shown with the same logarithmic 
stretch from contrasts of to 10^^. The single subtractions suppress the speckles outside 
the dark hole by a factor of 5 - 50. The bright point sources outside the dark hole are ghosts 
produced by aliasing outside the control radius of the DM. 
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Fig. 3. — Left:RAW, Right:NSDI. Gallery of single wavelength images and single differences 
with a nominal contrast level of 10~^. All images are shown with the same logarithmic 
stretch from contrasts of to 10^^. The single subtractions suppress the speckles outside 
the dark hole by a factor of 5 - 50. The bright point sources outside the dark hole are ghosts 
produced by aliasing outside the control radius of the DM. 
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Fig. 4. — Left:RAW, Right:NSDI. Gallery of single wavelength images and single differences 
with a nominal contrast level of 10~®. All images are shown with the same logarithmic 
stretch from contrasts of to 10^^. The single subtractions suppress the speckles outside 
the dark hole by a factor of 5 - 50. The bright point sources outside the dark hole are ghosts 
produced by aliasing outside the control radius of the DM. 



Fig. 5. — Left:RAW, Right:NSDI. Gallery of single wavelength images and single differences 
with a nominal contrast level of 10~^. All images are shown with the same logarithmic 
stretch from contrasts of to 10^^. The single subtractions suppress the speckles outside 
the dark hole by a factor of 5 - 50. The bright point sources outside the dark hole are ghosts 
produced by aliasing outside the control radius of the DM. 
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Fig. 6.— Double differenced images - (F4 - F3) - (F3 - F2) ((F3 - F2) - (F2 - Fl)) and 
(F5 - F3) - (F3 - Fl) for a nominal contrast level of 10~^. These images are shown with the 
same logarithmic stretch from contrasts of to 10~^ as Fig. [21 Significant speckle residuals 
are present compared to the single differenced images. In contrast to theory (Marois et al. 
2000), the double difference method is less effective at suppressing speckle noise than a single 
difference. The bright point sources outside the dark hole are ghosts produced by aliasing 
outside the control radius of the DM. 
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10"! A 




-40 -20 20 40 

AX (nm) from 798.9 nm 

Fig. 7. — Single image Contrast as a function of wavelength (AA in nm from 798.9 nm) 
inside the dark hole. Contrast is degraded by a factor of >2 at wavelength differences of 
±30 nm from 798.9 nm. The degradation of contrast as a function of wavelength is most 
pronounced at the highest nulling wavelength contrast level (10~^). 
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Fig. 8. — Single image Contrast as a function of wavelength (AA in nm from 798.9 nm) inside 
the dark hole (10~^ average contrast) and from model predictions. Contrast is degraded by a 
factor of >2 at wavelength differences of ±30 nm from 798.9 nm. The degradation of contrast 
as a function of wavelength is most pronounced at the highest nulling wavelength contrast 
level (10~^). Model predictions for this experiment produce a qualitatively similar result. 
These models are for slightly higher intrinsic contrast levels calculate contrast in somewhat 
different regions than used in the experiment, which is why there is some discrepancy between 
model and experiment results. The 4-5 A/D model samples only the deepest parts of the 
dark hole, while the 4-10 A/D model samples a much larger and somewhat less corrected 
region, hence, the observed degradation in contrast between the models. 



